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We consider a spin-orbit coupled configuration of spin- 1/2 interacting bosons with equal Rashba 
and Dresselhaus couplings. The phase diagram of the system is discussed with special emphasis 
to the role of the interaction treated in the mean- field approximation. For a critical value of the 
density and of the Raman coupling we predict the occurrence of a characteristic tri-critical point 
separating the spin mixed, the phase separated and the single minimum states of the Bose gas. 
The corresponding quantum phases are investigated analyzing the momentum distribution, the 
longitudinal and transverse spin-polarization and the emergence of density fringes. The effect of 
harmonic trapping as well as the role of the breaking of spin symmetry in the interaction Hamiltonian 
are also discussed. 
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A large number of papers has been recently devoted to 
the theoretical study of artificial gauge fields in ultracold 
atomic gases (for a recent review see, for example, 
First experimental realizations of these novel configura- 
tions have been already become available [2|, |3| . This field 
of research looks very promising from both the theoreti- 
cal and experimental point of view, due to the possibility 
of realizing exotic configurations of non trivial topology 
with the emergence of new quantum phases in both 
bosonic Q and fermionic gases, and the possibility to 
simulate electronic phenomena of solid state physics. In 
the case of Bose gases a key feature of these new systems 
is the possibility of revealing Bose-Einstein condensation 
in single-particle states with non-zero momentum. 

In the paper [3] it was possible to realize spin-orbit 
coupling in Bose-Einstein condensates, and to identify 
experimentally a quantum phase transition from a spa- 
tially spin mixed state to a phase-separated configuration 
by tuning the Raman coupling between two hyperfine 
states of 87 Rb atoms. Important features of the resulting 
phases were anticipated in the paper by Ho and Zhang 
[3] and discussed in detail in the same experimental pa- 
per 0. The purpose of this letter is to provide a the- 
oretical description of the phase diagram corresponding 
to the spin-orbit coupled Hamiltonian employed in jjjj. 
We point out the occurrence of an important density de- 
pendence in the phase diagram which shows up in the 
appearance of a tri-critical point that, to our knowledge, 
has never been predicted for such systems. 

We will consider the mean-field energy functional (for 
simplicity we set Ti = m = 1) 
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describing an interacting spin- 1/2 Bose-Einstein conden- 
sate, where ip a and ipb are the condensate wave functions 
relative to the two spin components interacting with the 



coupling constants gij = Aitdij , with the correspond- 
ing s-wave scattering lengths, and 



h = 



(p x -k a z ) + p 2 ± 
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is the single-particle Hamiltonian characterized by equal 
contributions of Rashba [8j and Dresselhaus [9] spin-orbit 
couplings and a uniform magnetic field in the x-z plane. 
In Eq.(j2]) Q is the Raman coupling constant accounting 
for the transition between the two spin states, fco is the 
strength associated with the spin-orbit coupling fixed by 
the momentum transfer of the two Raman lasers, S fixes 
the energy difference between the two single-particle spin 
states, <Ji are the usual 2x2 Pauli matrices, while V ex t 
is the external trapping potential. 

In the first part of the paper we will consider uniform 
configurations, neglecting the effect of the trapping po- 
tential (T4xt = 0) an d assume a spin symmetric interac- 
tion with g aa = gbb = g and 6 = 0. The ground state 
condensate wave function will be determined using a vari- 
ational procedure based on the following ansatz for the 
spinor wave function: 
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where TV is the total number of atoms, V is the vol- 
ume of the system. For a given value of the average 
density n = N/V, the variational parameters are then 
Ci, C2, k\ and 0. Their values are determined by mini- 
mizing the energy (pQ) with the normalization constraint 
Ei=„, b S<Pr\A\ 2 = N (i.e. |d| 2 + |C 2 | 2 = 1). Mini- 
mization with respect to yields the general relationship 
= arccos(/ci/&o)/2 (0 < < tt/4:), fixed by the single- 
particle Hamiltonian (j2j). Once the other variational pa- 
rameters are determined, one can calculate key physical 
quantities like, for example, the momentum distribution 
accounted for by the parameter fci, the longitudinal and 
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transverse spin polarization of the gas 



(a z ) = |(|C 1 | 2 -|C 2 | 2 ), (^) = -VM3 



and the density 
n{x) = n 



l + 2\C 1 C 2 \^\-^-cos(2k 1 x 



(4) 
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where <j) is the relative phase between C\ and G 2 . The 
ansatz (J3j) describes exactly the ground state of the 
single-particle Hamiltonian ho (ideal Bose gas). In this 
case, for tt < 2fcg, the energy, as a function of k\ , exhibits 
two minima located at the values ±fco\A — ^ 2 /4/cq and 
the ground state is degenerate, the energy being indepen- 
dent of the actual values of C\ and C 2 . For Q > 2/cq the 
two minima disappear and all the atoms condense into 
the zero momentum state k\ = 0. 

The same ansatz is well suited to discuss the role of 
interactions. By inserting (|3]) into (pQ), we find that the 
energy per particle e = E/N takes the form 
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y /l$-k*-F(l3)^+G 1 (l + 2(J) (6) 



where we have defined the dimesnionless parameter j3 
\Ci \ 2 \C 2 \ 2 (0 < P < 1/4), and the function 



F(fi) = (k* - 2G 2 ) + 4 (Gi + 2G 2 ) P 



(7) 



with the interaction parameters G\ — n (g + g a b) /4, 
G2 = n (g — g a b) /4. The variational parameters to min- 
imize the energy are then k\ and /3. 

Let us first consider minimization with respect to fci. 
If Q > 2F(j3>) the energy ((6]) is an increasing function of 
k\ and the minimum takes place at k\ = 0. If instead 
£7 < 2F(/3) one finds that s is minimized by the choice 



ki(P) = k 0] 



n 2 



4[F(/3)] 
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which generalizes the ideal gas result (F(/3) = k%). 
Eqs. ((TJ)-(IHJ) explicitly show that the momentum distribu- 
tion is modified by the interactions. We find the following 
result for the energy per particle 



8F(fi) 



Gi + G 2 (1 - 4/3) . 
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The ground state of the system can be found by looking 
for the minimum of (|9j) with respect to j3. One can easily 
prove that the second order derivative of ((9]) with respect 
to j3 is negative. The ground state is then compatible with 
the three following phases: 

(I) The spin mixed or "stripe" phase with k\ ^ 0, /? = 
1/4 and hence (<j z ) = 0. In this phase the atoms condense 
in a superposition of two plane wave states with wave 



vector ±ki and the density (j5j) exhibits fringes. This 
configuration is characterized by a degeneracy associated 
with the relative phase between the coefficients C\ and 
C2 which fixes the actual spatial position of stripes. 

(II) The separated phase with k\ 7^ 0, j3 — and hence 
(<j z ) 7^ 0, where the atoms condense into a single plane 
wave state with wave vector either k\ (C2 = 0) or — k\ 
(Ci = 0), the actual value being determined by a mech- 
anism of spontaneous spin symmetry breaking. 

(III) The single minimum or "zero momentum" phase 
with k\ = and (a z ) = where the atoms condense in 
the zero momentum state. 

We first notice that the spin mixed phase is compatible 
only with positive values of the interaction parameter 
G2, favoring anti- ferromagnetic configurations. In fact 
in the opposite case G2 < 0, the first order derivative 
de/dj3 is always positive and the ground state is always 
in the phase separated configuration (II) or in the "zero 
momentum" phase (III). 

In the most interesting G2 > case, the system will 
be always in the phase (I) for small values of the Raman 
coupling constant Q. If the condition 



kl > 4G 2 



4Gf 



(10) 



is satisfied, the systems will exhibit a phase transition (I) 
to (II) at the frequency 
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This generalizes the result derived in [7[, which corre- 
sponds to the low density (or weak coupling) limit of 
(fTTj) . i.e. Gi, G2 <C &q. The transition frequency in this 
limit approaches the density independent value 



ft (I - n) - 



2fcgv/2 7 /(l + 2 7 ) 
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where we have introduced the dimensionless interaction 
parameter 7 = G 2 /Gi = (g - g a b)/(g + 9 ah)- By further 
increasing ^, the system will enter the phase (III) at the 
frequency 



n< n - in > = 2 {kl - 2G 2 ) 



(13) 



If instead the condition (flQ|) is not satisfied, the transi- 
tion will occur directly from the phase (I) to (III) at the 
frequency 



l^ 1 - 1 ") = 2 {kl + G 1 ) - 2 [(fc 2 + Gi) Gi] 



1/2 
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In the strong coupling limit G\ ^> k^ ([T4j) approaches 



the constant value k^. 



The critical point where the phase (II) disappears is 
fixed (see Eq. ffTQl) ) by the condition G{ c) = fcg/47 (1 + 7), 
corresponding to the critical value 



n^ = kl/{2 lg ) 



(15) 
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FIG. 1: (Color online) fci, energy per particle E/N, transverse 
and longitudinal spin polarization (<j x ) and |{<7 Z )| as a func- 
tion of Q. Red solid lines: "stripe phase" k\ ^ and f3 — 1/4; 
blue solid lines: "separated phase" ki ^ and f3 — 0; green 
solid lines: "zero momentum phase" k\ — 0; open circles: 
ground state. The parameters: Gi/ko = 0.2, G^/fco — 0.05 
(a)-(d), G 2 /fc 2 =0.16 (e)-(h). 



FIG. 2: (Color online) Spin polarization \(cr z )\ (a) and ki/ko 
(b) as a function of Q and density n/n^ in three different 
phases. The white solid lines represent the phase transition 
(I II), (II-III) and (I-III). The parameters: g = 100 a s , where 
as is Bohr radius, 7 = 0.0012, k% = 27rx80Hz, corresponding 
to the critical density n (c) = 4.37 x 10 15 cm -3 . 



for the density. If n < n^ c \ one has two transitions (I-II 
and II-III), while if n > n^ c \ only one phase transition 
(I-III) can take place. 

In FigJTJ we plot the momentum fci, the energy per 
particle E/N, the transverse and longitudinal spin po- 
larizations (a x ) and \(cr z )\ as a function of Q for n < 
(left column) and for n > (right column). In ad- 
dition to the results for the ground state (open circles), 
we also show the various quantities for the three phases 
(colorful lines). Figure (a)-(d) reveals the emergence of 
the phase transitions (I-II) and (II-III), while in (e)-(h) 
there is only the transition (I-III). The figures also show 
that the transitions (I-II) and (I-III) are accompanied by 
a jump in fci (see (a) and (e)) and consequently in (a x ) 
(see (c) and (g)). In particular the jump in k\ associated 
with the transition (I-III) is sizable and should be easily 
observable in experiments. On the other hand only the 
transition (I-II) is accompanied by a jump in the longi- 
tudinal spin polarization \ (cr z )\. The transition (II-III) is 
instead characterized by a continuous behavior of the rel- 
evant physical parameters. The experimental conditions 
of [3[ correspond to values of the average density n much 
smaller than so that the jump in fci could not be 
detected being too small at the transition (I-II). On the 
other hand the occurrence of this phase transition was 
clearly revealed by the analysis of the spin distribution 
after time of flight (see Fig.2c of [3[). 

In Fig|2]we show the phase diagram for three differen 



phases. The value of the spin polarization \{cr z )\ an d fci 
are reported in (a) and (b) respectively. The transition 
lines separating different phases merge at a tri-critical 
point at n = n^ c \ The value of \(cr z )\ always vanishes 
for n > n^ c \ However the phase transition (I-III) is well 
identified by the behavior of the momentum k± . The pa- 
rameters employed in Figj2] correspond to rather large 
values of the critical density. More accessible values of 
can be obtained employing smaller values of fco. Ex- 
perimentally this can be achieved by playing with the an- 
gle between the two Raman lasers. Reducing the value 
of fco nas a l so the advantage of increasing the spacial 
separation between the fringes in the "stripe phase" (I), 
thereby making their experimental detection easier. 

Let us now discuss the effect of the trap. In order to 
simplify the analysis we have considered harmonic trap- 
ping with frequency uoq only along the x-axis. Without 
interaction, one can calculate the ground state using a 
similar variation ansatz, replacing the plane waves in (j3j) 
by the functions e ±lklX e~^ x / 2 , corresponding, in the 
absence of the gauge field, to the usual harmonic oscil- 
lator Gaussians. The energy per particle is easily calcu- 
lated and reads: 

up fcg - k\ Q r~ 2 — 

p ( 16 ) 

- (C*C 2 + C*d) ft e-^ 
2fc 
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FIG. 3: (Color online) Spin polarization \{cr z )\ as a function 
of Q for the trapped case (red solid line) , and for the uniform 
case using the density in the center of the trap (blue dashed 
line). The parameters are chosen as follows: cjo = 27T x 20 Hz, 
ko/uoo — 4, g aa — gbb — 101.20 as, g a b — 100.99 as, where as 
is Bohr radius. The density in the center of trap corresponds 
ton - 3.9 x 10 13 cm" 3 . 

The ground state can be found by minimizing e with re- 
spect to fci, C\ and C2 with the normalization constraint. 
The first term in ([T6]) is just the zero point energy due 
to the presence of the trap. The following two terms are 
the same as for the uniform case without interactions, i.e. 
(|6j) with G\ = G2 = 0. The last term shows the effect 
of the trap, fixing the relative phase between the coeffi- 
cients C\ and C2 in the ground state. Consequently the 
degeneracy occurring in the uniform case will be lifted 
even in the absence of interactions (where = 0). Phys- 
ically this is the consequence of the non orthogonality of 
the two Gaussians states centered in dzfei. According to 
([T6]h for ki ^ 0, the system prefers to stay in the spin 
mixed phase, and exhibits density modulation in space 
even without interactions. On the other hand the interac- 
tion is crucial for the appearance of the phase separated 
configuration. Since the last term of (fT6|) scales expo- 
nentially, the effect of the trap is weak for k\ ^> co>o, and 
becomes more and more important when k\ is compara- 
ble to UOq. 

To describe the role of the interaction we implement 
the mean-field approximation by solving numerically the 
Gross-Pitaevskii equation for the condensate wave func- 
tion using the gradient method in the same ID trapping 
conditions. We find that the properties discussed in the 
first part of the work for the uniform system almost hold 
in the trapped case. In Fig|3]we show an example of the 
numerical calculation. The spin polarization in presence 
of trap as a function of Q (red solid line) is compared with 
our analytical results for the uniform case (blue dashed 
line), using the density in the center of the trap. There 
is good agreement between the two curves. 

We finally discuss the case 5 =^ and g aa 7^ #55, cor- 
responding to broken spin symmetry. In general for uni- 
form case, one can introduce three interaction parame- 
ters: Gi = n(g aa + g bb + 2g ab )/8, G 2 = n(g aa + g bb - 
2g a b)/8, and G 3 = n(g aa - gbb)/^- In the case of 87 Rb 
atoms, the scattering lengths relative to the spin states 



\F = 1, nip = 0) and \F = 1, nip = —1) are usually 
parameterized as a aa = Co, a^b = Co + C2 = a a 5, with 
c = 7.79 x 10- 12 Hz cm 3 and c 2 = -3.61 x lO" 14 Hz cm 3 . 
This corresponds to G2 = G3 <C G\. However, since the 
differences among the scattering lengths are very small, 
by properly choosing the detuning <$, this effect can be 
well compensated, and the properties of the ground state 
remain the same as before. For example, using first order 
perturbation theory, one finds 

e (i) = ^ 3 + 0|(|C 1 | 2 -|C 2 | 2 ) (17) 

where we have considered the low density (weak cou- 
pling) limit. By choosing S = — 2G3 the correction (fT7|) 
identically vanishes and the transition frequency (fT2j) is 
not consequently affected by the inclusion of the new 
terms in the Hamiltonian. Using the 8T Rb parameters 
introduced above we then find the value f^ D n) = 0.19 E L 
(El = in agreement with the findings of [3] corre- 

sponding to n/n^ <C 1. For higher densities, the value 
of S should be instead taken ft dependent in order to 
ensure exact compensation. 

In conclusion, we have investigated the phase diagram 
of a spin-orbit coupled two-component Bose-Einstein 
condensates using a variation ansatz based on the mean- 
field approximation. We predict a rich phase diagram 
characterized by the occurrence of three different quan- 
tum phases, and by a characteristic tri-critical point 
where the three phases merge at a critical value of the 
density and of the Raman frequency. 
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